Azimuthal anisotropy: transition from hydrodynamic flow to jet suppression 
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Measured 2 nd and 4 th azimuthal anisotropy coefficients V2,4,(N pa ,rt,PT) are scaled with the initial 
eccentricity £2,4(A r part ) of the collision zone and studied as a function of the number of participants 
-/Vpart and the transverse momenta pr- Scaling violations are observed for pr < 3 GeV/c, consistent 
with a p\ dependence of viscous corrections and a linear increase of the relaxation time with pr- 
These empirical viscous corrections to flow and the thermal distribution function at freeze-out 
constrain estimates of the specific viscosity and the freeze-out temperature for two different models 
for the initial collision geometry. The apparent viscous corrections exhibit a sharp maximum for 
Pt > 3 GeV/c, suggesting a breakdown of the hydrodynamic ansatz and the onset of a change from 
flow-driven to suppression-driven anisotropy. 
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A central objective of the experimental program at the 
Relativistic Heavy Ion Collider (RHIC) is to delineate 
the thermodynamic and transport properties of the hot 
and dense matter produced in energetic heavy ion colli- 
sions. This matter can equilibrate to form a hot plasma 
of quarks and gluons (QGP) [[J which rapidly expands, 
cools, and then hadronizc to produce the observed par- 
ticles. The hydrodynamic-likc expansion of the QGP, as 
well as its interactions with hard scattered partons, re- 
sults in the anisotropic emission of hadrons relative to 
the reaction plane. At mid-rapidity, the magnitude of 
this momentum anisotropy is characterized by the even- 
order Fourier coefficients; v n = (V n ( A( ^) , n = 2,4,..., 
where A(f> is the azimuth of an emitted hadron about 
the reaction plane, and brackets denote averaging over 
particles and events. 

The coefficients for hadrons with low transverse mo- 
menta (pr < 2 GeV/c) can be understood in terms of 
flow or partonic interactions which drive pressure gra- 
dients in an initial "almond-shaped" collision zone pro- 
duced in non-central collisions |i Q. For higher trans- 
verse momenta (j>t > 5 GeV/c) the coefficients can be 
attributed to jet quenching [ - the process by which 
hard scattered partons interact and lose energy in the hot 
and dense QGP, prior to fragmenting into hadrons. This 
energy loss manifests as a suppression of hadron yields 
|J which depends on the average distance that partons 
propagate through the QGP. Thus, «2,4 stem from the 
fact that partons which traverse the QGP medium in a 
direction parallel (perpendicular) to the reaction plane 
result in less (more) suppression due to shorter (longer) 



parton propagation lengths [| J . This path- length de- 
pendence is exemplified in the recently observed scaling 
patterns for hadron suppression pAL O. The transition 
from flow-driven to suppression-driven anisotropy is still 
poorly understood, and it remains a challenge to find a 
single consistent theoretical framework that gives an ex- 
planation of t>2,4 measurements over the full pt range. 

The magnitude of i>2,4, as well as their detailed de- 
pendence on pt and collision centrality (or number of 
participants iVp art ), give invaluable insights on the ther- 
modynamic and transport coefficients of the QGP. In par- 
ticular, flow measurements (i>2(pt) and W2(-/V part )) have 
been used to estimate the specific shear viscosity (i.e., 
the ratio of shear viscosity rj to entropy density s of the 
plasma) via comparisons to viscous relativistic hydrody- 
namic calculations p jj] . The reliability of these ^ es- 
timates is influenced not only by the uncertainties in the 
initial conditions for hydrodynamic evolution, but also 
by ambiguities in the departure from equilibrium on the 
freeze-out surface. For a viscous fluid, this distribution 
(/) is of the form 



dN 



dyprdpTi 



/o + Sf = /o 1 + G 



Tf 



(1) 



where /o is the equilibrium distribution, Tf is the freeze- 
out temperature, C ~ „ n T and a ranges between and 



69 1 llH, pi ; r is the time scale of the expansion. The factor 

70 (5/, which results from a finite shear viscosity, is known to 

71 dominate the calculated viscous corrections to V2{pt) for 

72 Pt ^ 1 GeV/c (Hj. However, its momentum dependence 
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to extract estimates 

of the viscous corrections to V2,i{px, ^Vpart)- I n turn, we 

95 use these estimates to explore the px dependence of Sf 

96 and the transition from flow-driven to suppression-driven 

97 anisotropy. We find viscous correction factors for pt < 3 

98 GeV/c that validate the commonly assumed p T depen- 

99 dence of Sf, and give constrained estimates for - and 
ioo the freeze-out temperature 2). For px > 3 GeV/c, the 
ioi apparent viscous corrections signal a breakdown of the 

102 hydrodynamic ansatz. 

103 The 7J2.4 data employed in our analysis are selected 

104 from the high-precision PHENIX measurements recently 
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FIG. 1. (color online) Comparison of ^2/22 vs. -/V par t for sev- 
eral pr selections, obtained from perfect fluid (a) and viscous 
hydrodynamic (b) simulations ||19| of Au+Au collisions. 



73 and associated relaxation time tr{jp) is not known a pri- 

74 or% and it is unclear whether or not it is proportional to 

75 p T {a — and tr oc p) as has been generally assumed 

76 in hydrodynamic calculations 1 17, IfJ. The freeze-out 

77 temperature and the px for which large viscous correc- 

78 tions render a breakdown of viscous hydrodynamics are 

79 also not well established experimentally. 

so The influence of viscous corrections on the eccentricity- 
8i scaled anisotropy coefficient ^ (tv'*' ) i s illustrated in 

82 Figs. (a) and (b) where the results of hydrodynamic 

83 simulations (with the code of Dusling and Teaney J 5 _) ) 

84 are shown for | = and 0.2, respectively. Fig. [I (b) 

85 shows that viscous effects reduce V2(N part ,px) and break 
se the scale invariance of ideal hydrodynamics evidenced in 

87 Fig. [P (a), i.e. there are deviations away from the essen- 

88 tially flat N palt dependence expected for ideal hydrody- 

89 namic scaling. These deviations from eccentricity-scaling 

90 can be used to estimate and characterize viscous correc- 

91 tions |[g 2|. 

92 Here, we use the eccentricity scaled anisotropy coeffi- 

«2(A f part,PT) J «4 (Npart ,PT ) 



2. (color online) Comparison of ^2/22 vs 
54 vs 



A^part (a) and 
-/Vpart (b) for several pr selections as indicated pj. 



FIG 

114/. 

The dashed curves indicate a simultaneous fit to the data in 
(a) and (b) [for each px] with Eq. (y). The V2 4 data are from 
Ref. @. 



105 reported for Au+Au collisions at y / 5jv"/v = 200 GeV [Q. 
we These data show that both V2 and V4 have a strong de- 

107 pendence on px and centrality. The large increase in 

108 V2,i(px) from central to peripheral events is especially 
no important to our study of viscous corrections. 

in Monte Carlo (MC) simulations |2 j| were used to calcu- 

112 late the event-averaged geometric quantities used in our 

113 analysis. For each collision, the values for 7V part and the 

114 number of binary collisions iV C oii were determined via the 

115 Glauber ansatz [||(J . The associated values for the trans- 
it verse size R, area S and eccentricities £2,4 were then cval- 
117 uatcd from the two-dimensional density of sources in the 
us transverse plane /O s (rj_) using two principal models; a 

119 modified version of the MC-Glauber approach | | and 

120 the factorized Kharzeev-Levin-Nardi (MC-KLN) model 

12, nfn 

122 For each event, we compute an event-shape vector S n 

123 and the azimuth of rotation VP* for the n-th harmonic of 

124 the shape profile |]j; 

S nx = S n cos(n^) = / drj_p s (rjJw(rjJcos(nr/>), 



S ny = 5„sin(n**) = / dr±p s (r±)uj(r±)sm(n<j)), 



11 \ >->nx 



where (f> is the azimuthal angle of each source and the 
weight oj(t±) = r± 2 ; £2,4 were calculated as: 

£ 2 = (cos 2(0 -*5)), £4 = (cos A{<j> -*5)>, (2) 
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FIG. 3. (color online) v^lei vs. iVpart for several (pr) values 
as indicated. The filled circles are the same as in Fig. ua). 
The dashed curves show fits to the data obtained with Eq. 
(|). The v 2 data are from Ref. (Jl. 



127 where the brackets denote averaging over sources, as well 

128 as events belonging to a particular centrality or impact 

129 parameter range. For the MC-Glauber calculations, an 

130 additional entropy-density weight was applied reflecting 

131 the combination of spatial coordinates of participating 

132 nucleons and binary collisions l |; 
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FIG. 4. (color online) K* vs. (p T ) ( a ) an d (b), and KR vs. 
Apart (c) and (d), extracted with MC-KLN (left panels) and 
MC-Glauber geometry (right panels). The filled circles in (a) 
indicate results from the simultaneous fits shown in Fig. 0. 
The dashed curves in (a) and (b) show a fit to the data for 
{pt) 2 < 10 [GeV/c] 2 . 



158 In lieu of detailed model comparisons p2| , we estimate 

159 the magnitude of the viscous corrections by parametriz- 
leo ing the observed scaling violations with a Knudsen num- 
i6i ber (K = X/R) ansatz akin to that in Refs. (20, 21 ; 



133 where a± = 0.14 was constrained by multiplicity mea- 

134 surements as a function of IVpart f° r Au+Au collisions 

135 |J0J. Note that £2,4 (cf. Eq. (0)) correspond to ^2,4 mea- 

136 surements in the so-called participant plane fj; this is 

137 analogous to the measurement of ^2,4 with respect to the 
OS 2 nd order event-plane as described in Ref. 1 3] . A corre- 
139 lation between the principal axes of the quadrupole (^2) 
wo and hexadecapole (^|) density profiles was also intro- 

141 duced to account for contributions to 7J4 from v^ \~\- 

142 This correlation has a significant influence only on £4 in 

143 the most central collisions |_) . 

144 Figures g and y show eccentricity-scaled values of 

145 V2,a(pt, A^part) obtained with MC-KLN eccentricities for 

146 several px cuts. The low-px selections are almost flat, 

147 i.e. small scaling violations. These violations gradu- 

148 ally increase with pt over the pt range indicated in Fig. 

149 g. That is, the data points slope upward progressively 

150 (from low to high -/V pa rt) as the (pr) is increased. Fig. 

151 U shows that this trend reverses to give scaling viola- 

152 tions which decrease with increasing (pr), for (pr) > 3 

153 GeV/c. This inversion could be an indication for the on- 

154 set of suppression-driven anisotropy as discussed below. 

155 Similar scaling performed with MC-Glauber eccentrici- 

156 ties, show the same trends exhibited in Figs. and w, 

157 albeit with larger scaling violations, as discussed below. 



V2k(Pr) _ V2 k (Pr) 



£2fc 



2k\ 
£2k 



1 



1 '■■ 



K*(pt) 



Ka 



fc = l,2, 



(3) 



where K*(px) characterizes the magnitude of the viscous 



correction for a given pt, 



i(Pr) 



are the eccentricity- 



164 scaled coefficients expected from ideal hydrodynamics, A 

165 is the mean-free path, and K is a constant estimated to 

166 be 0.7 ± 0.03 with the aid of a transport model 

167 For each pr selection, [K*(px)]~ 1 



( dN 
\ dy 



P<PT)$%i 



Npart) is evaluated by fitting 



T>2,4(PT, -Npart) 
!,4(JVpart) 



VS. 

Ap^t with Eq. (§) (cf. curves in Figs |"and'|j! The 

170 fit parameters (3(pt), so obtained, allow the determina- 

171 tion of K*(pt) as a function of A par t. Note that a model 

172 uncertainty in the value of Kq would lead to an accompa- 

173 nying uncertainty in the magnitude of K* . However, the 

174 consistency of our procedure with hydrodynamic models 

175 has been tested via fits to P ?* T P T vs. N nart . obtained 

176 for specified values of 2 [F^j. These fits lead to ^ values 

177 which reproduce the input values to the hydrodynamic 

178 simulations. 
Figures (a) and (b) show the values of K* vs. p^, ex- 



tracted for N, 



351 with MC-KLN and MC-Glauber 
geometries respectively (the plots for other values of 



part 



182 iVp ar t show similar trends, but with different intercepts). 

183 The filled circles in Fig. [J (a) show results from the simul- 

184 taneous fits indicated by the dashed curves in Figs. || (a) 

185 and (b). The squares show results for fits which employ 

186 only " 2 e 7^ pa \ data J2JJ; both are in good agreement. 

187 The K* values shown in Figs. H (a) and (b) indicate a 

188 linear dependence on p T , for (pr) < 10 [GeV/c] 2 , that 

189 demonstrates a non-zero viscosity and the p T dependence 
iso of 5f commonly assumed in hydrodynamic simulations 
i9i p l]], |jj|. In contrast, the data for 10 < (p T ) 2 < 40 

192 [GeV/c] 2 show a striking trend inversion. We interpret 

193 this as a signal for the breakdown of the hydrodynamic 

194 ansatz when K* ~ 1, as well as an indication for the 

195 onset of suppression-driven anisotropy. Note that such 

196 a scenario would lead to improved eccentricity scaling 

197 for pt > 3 GeV/c because the eccentricity encodes the 

198 variation of the path length relative to the orientation of 

199 the reaction plane, and the suppression of hadron yields 

200 has been found to increase as -h= for a similar pt range 

201 ||. 

202 To obtain estimates for K and 2} for each value of 



P 2 i 



-Npart, we use our observation that K*(pt) 
conjunction with the first order expansion of ^(pt 

f*^ dA(t> cos(2A0) 



dy pf dpf dA0 
H 3 M 



to obtain 



2„5 <cos(2A$) PT = 

•I— 7r " dy pi dpf dill) 

206 the expression K*(pt) = K + ^- (^r) ; the constant 

207 B was cross-checked via fits to the results from hydro- 

208 dynamic simulations. Fits to K* vs. p T were performed 

209 with this fit function. The dashed curve in Fig. 0(a) indi- 

210 cates such a fit for (px) < 10 [GeV/c] 2 ; it gives the val- 
2u ues K = 0.09±0.01 (from the intercept) and T f = 162±11 

212 MeV (from the slope). The same fit to the K* values in 

213 Fig. [J(b) (extracted with MC-Glauber eccentricities), 

214 give the values K = 0.17 ± 0.007 and T f = 173 ± 11 



235 Estimates for - were obtained via the expression 

236 2 m XTc s = (RKTc s ) where the sound speed c s = 

237 0.47±0.03 c was obtained from lattice calculations p J for 

238 the mean temperature T = 220 ±20 MeV Q. This gives 

239 the estimates 4tt2 = 1.1 ± 0.1 and 4tt2 = 2.1 ± 0.2 for 

s s 

240 the K values extracted using MC-KLN and MC-Glauber 

241 eccentricities [respectively] in central and mid-central col- 

242 lisions. These estimates are in agreement with the low 

243 value from prior extractions |i 17], |l9|, ^l] ^3|, || ] . 

244 In summary, we have used eccentricity scaled 

245 anisotropy coefficients to extract estimates of the 

246 strength and role of the viscous corrections. These es- 

247 timates show a quadratic increase with pt (for px < 3 

248 GeV/c) that validates a non-zero viscosity and a relax- 

249 ation time which grows with pt- The extracted vis- 

250 cous corrections also constrain the estimates 4-7T- ~ 
25i l.liO.l (2.1±0.2) and T f = 162±11 MeV (173±11 MeV) 

252 for MC-KLN (MC-Glauber) collision geometries for a 

253 strongly coupled plasma. The onset of a transition from 

254 flow-driven to suppression-driven anisotropy is signaled 

255 by a sharp maximum of the apparent viscous corrections 

256 for pt > 3 GeV/c. These results provide valuable con- 

257 straints for input parameters to more detailed viscous 

258 hydrodynamic calculations. 
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